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Gribov's theory of light quark confinement implies the existence of two kinds of 
scalar bound states. The phase diagram of the three-flavor QCD is mapped out 
in the (m^ — m^)— plane with help of the SUl{3) x SUr(3) linear sigma model 
supplemented with the assumption that the masses of the so-called superbound 
scalars do not change under the variation of the pion and kaon mass. The phase 
boundary along the m n = mx line is found in the interval 15 MeV < m cr ; t < 
25 MeV, irrespective which fo — a linear combination is identified with the pure 
superbound state. 



1. Introduction 

Volodia Gribov repeatedly has expressed in the early nineteen-eighties his 
view that deconfinement in QCD is going to turn out to be a complete 
analogue of atomic ionisation, which one would not call a phase transition. 
There is increasing evidence that the temperature driven transition from 
the hadronic phase to quark-gluon plasma is a smooth crossover indeed . 
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He based his confinement theory also on an analogy with a phenomenon 
of atomic physics, e.g. supercritical binding. This would occur for light 
quarks irrespective to the accurate value of their mass 2 . The role of chiral 
symmetry is not evidently important in this theory, and the consequences 
on its spontaneous breakdown were not yet fully clarified 3 . 

In this contribution to the commemoration of the 75th anniversary of 
V.N. Gribov, we exploit semi-quantitatively his ideas on the nature of the 
scalar meson sector when discussing the chiral symmetry restoring quark- 
hadron transition in different points of the — m^-plane with help of 
a variant of the linear sigma model {LaM). In a recent publication 4 we 
constructed a continuation of the parameters of this model from the physical 
point to an arbitrary mass-point m^^niK requiring agreement with results 
of Chiral Perturbation Theory (ChPT) for the decay constants f n , fx and 
the trace of the squared mass matrix in the 77 — 77' sector, M 2 , at the 
tree level. In the 77 — 77' sector the predictions of the two models for the 
separate mass eigenvalues nicely coincide. It was emphasised that for a 
complete specification of the parameters of the model one needs a single 
extra information on the dependence of scalar spectra on the masses of the 
pseudoscalar nonet. This information is beyond ChPT, therefore in general 
only ad hoc assumptions can be made and tested through the consequences. 

According to Gribov's confinement theory in addition to normal qq 
states also superbound states do exist, which contain quarks with negative 
kinetic energy. Some repulsive interaction between such quarks increases 
their energy above zero and makes these associations physically observ- 
able 5 . These states are of smaller size and higher energy than the normal 
scalar meson states, therefore it was suggested to associate them with fo in 
the isoscalar and with a in the isovector channel. In the SU{2) x SU{2) 
model they form another 0(4) quartet in addition to <r, ir. It is rather nat- 
ural to assume that the masses of the superbound states are not sensitive 
to the variation of m^. 

In the three flavor case a similar "doubling" of the multiplct structure 
can be assumed, if the strange quark is light enough. In this paper we shall 
explore the boundary of the first order transition region around the chiral 
point in the three-flavor QCD 6 taking into account the extra requirement 
to have scalars in the spectra whose mass does not vary with m n ,mK- 
The consequences of identifying the physically observable scalars with some 
mixture of the pure normal and pure superbound states will be investigated. 

Numerical investigations of the chiral symmetry restoration were done in 
the framework of lattice QCD and systematically improved for the 3-flavor 
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degenerate case m u = m,d = m s ^ 0. The initial estimate for the critical 
pseudoscalar meson mass, m cr i t (diag) w 290 MeV 7 was seen to be reduced 
to 60 — 70 MeV 8 or may be to even further down 9 when finer lattices and 
improved lattice actions are used. Very recently de Forcrand and Philipscn 
reported an estimate m cr i t ~ (0.1 — 0.2)T C , which very conservatively means 
merit ~ 15 — 30MeV 10 . One should be conscious of the fact, that it is 
extremely difficult to reach continuum results in this mass regime. 

Effective models (linear or non-linear sigma models, Nambu-Jona- 
Lasinio model) represent another, in a sense complementary, approach to 
the study of the phase structure, which one expects to work the better 
the lighter quark masses are used n > 12 > 13 . It is surprising that only mod- 
erate effort was invested to date to improve the pioneering studies of the 
SU(3) x SU(3) linear sigma model by Meyer-Ortmanns and Schaefer 11 
which used a saddle point approximation valid in the limit of infinite num- 
ber of flavors, and derived m cr it (diag) < 51 MeV. An extension of their work 
to unequal pion and kaon masses was achieved by C. Schmidt . He found 
™crit(diag) = 47 MeV and a phase boundary approaching the m^-axis 
rather sharply. The phase boundary was calculated also by Lenaghan 15 
using the Hartrcc-approximation to the effective potential derived in CJT- 
formalism. For the complete determination of the couplings of the three- 
flavor chiral meson model he fixed the T = mass of the a particle in 
addition to the experimental mass spectra of the pseudoscalar sector. The 
emerging phase boundary is rather sensitive to this mass. The estimate for 
77i cr it(diag) which one can extract from Fig. 3 of 15 for m a = 900 MeV is 
compatible with n > 14 . 

Our method of parametrisation and solution of LcrM was described in 
detail in 4 . Therefore here we shall only review the set of the equations to 
be solved for the determination of the transition point. The changes arising 
from the implementation of the insensitivity of Gribov's superbound scalars 
to the variation of the fundamental masses will be emphasised. We conclude 
by giving the most characteristic features of the phase boundary. 



2. LcrM parametrisation consistent with ChPT 

We outline first, how one obtains with ChPT the m w , m^-dependence of 
the masses and decay constants in the pseudoscalar sector. The dependence 
of the pion and kaon masses as well as of their decay constants on the quark 
masses were determined in 16 for the SU (3) x SU(3) nonlinear sigma model: 
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Here A, q are related to the quark masses, / is the coupling of the non- 
linear sigma model, fips defines the so-called chiral logarithm for each 
pseudoscalar meson (PS) proportional to ln(mps/M ) with M = 4-7r/. 
The determination of the low energy chiral constants Li is discussed in 
depth in 4 . 

One inverts the first two equations with 0{f~ 2 ) accuracy and finds the 
following m w , mx -dependence for the decay constants: 



U = f 

Ik = f 



1 - -jj (2/^ + n K - Avn 2 {L 4 + L 5 ) - 8m 2 K L 4 ) 
1 /3 

1 ~ 72 ( 4^ + Mt? + 2mk) - 4m^L 4 - 4m 2 K (L 5 + 2L 4 ) 



(5) 



The extension to the U(3) x 17(3) ChPT is somewhat more compli- 
cated. It was worked out in 17 > 18 > 19 > 20 an(1 allows the determination of 
m I) (m T , ttik), rn v i (m n ,m,K)- For the parametrisation of LaM one more in- 
dependent relation can be obtained from the mixing ij — rj' sector, for which 
in 4 we have chosen the trace of the 2x2 squared mass matrix, denoted by 



m~ + to~, : 



Ml = 2m 2 K - + 2(2m 2 K + ml)(3v™ - v™) 



(6) 
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Also the constants were determined in 4 from the T = properties of 
the Goldstone particles. 

Next, we turn to the problem of parameterising LaM to provide spectra 
the closest possible to ChPT. The Lagrangian of the SU L (3) x SU R (3) 
symmetric linear sigma model with explicit symmetry breaking terms is 
given 21 by 

L(M) = itr^M^M + &M*M) - f l (tr(M^M)) 2 (7) 
- f 2 ti{M^M) 2 - g (det(M) + det(Aft)) 
where M is a complex 3x3 matrix, defined by the <7j scalar and Hi pseu- 



doscalar fields M 



1 

:= —= y^((7i+i7ri)A,, with \ : i = 1 . . . 8 the Gell-Mann 
V ^ 



i=0 



matrices and Ao := y |l. The last two terms of (8) break the symmetry 
explicitly, the possible isospin breaking term 63173 is not considered. 

A detailed analysis of the symmetry breaking patterns which might 
occur in the system described by this Lagrangian can be found in 12 . The 
field expectation values (a ),(a s ) both contain strange (= y) and non- 
strange (= x) components: 

x = (a/2(gt ) + (<78»/>/3, y = «<r ) - \/2(cx 8 »/V3. (8) 

With help of the PCAC relations and the tree level mass formulae (see 
Table 1, where the x — y basis is used instead of — 8) the following ex- 
pressions can be derived for the couplings of LaM: 



h 



fJ 



U, y = (2f K -U)/V2, (9) 
(6/g - 3U)m 2 K - (2f K + f^ml - 2(f K - U)M 2 

KfK-u)(&f 2 K -&.f K u+m) 

2f K m 2 K + 2(f K - f,)ml - (2f K - U)M 2 
V2(8/i-8/K/. + 3/ 2 ) 
M 2 ^-^ 2 +Ah{x 2 +y 2 ) 

= hi 2 + X( 2 / x - U) 2M(.U f K f + f K \. 

The sources e x = (V2e + eg)/-^, e y = (eo — V2es)/VS, which explicitly 
break chiral symmetry are determined with help of the Gell-Mann-Oakes- 
Renner relations: 

e x = mlx, e y = ^(m 2 K -ml)x + m 2 K y. (10) 
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Table 1. Squared masses of the pseudoscalar boson nonet and 
their parity partners 



ml 


= - M g + 2(2/i - 
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f / 2 )z 2 


f 4/i j/ 2 -2gy 
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= -Mg + 4/ix 2 


f 4(/i 4 
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Vxy 








XX 


= -M( 2 , + 6(2/! - 


f / 2 )* 2 


+ 4f iy 2 + 2gy 


mi 

°yy 


= -^+4/1^ 


+ 12(/l 


+ h)y 2 


mi. 

&xy 


= Sfixy + 232: 







Note: The expressions of the squared masses of parity partners 
having the same isospin and hypercharge appear in one block. 
Different isomultiplcts arc separated by double lines. In the lowest 
big block the matrix elements of the mixing in the 17 — 77' and cr — fo 
sectors are given in the x — y base. 



The requirement of the agreement of LaM with ChPT is fulfilled when 
the {m-n, m7f)-dependence of the couplings x, y, /2, g, M 2 is determined with 
help of Eqs.(5,6). Before the procedure just described was first proposed in 
Rcf. 4 m^—rriK mass tuning was taken into account only in (10). The quality 
of this not fully complete paramctrisation (e.g. only the combination M 2 
of /1 and /j,q is determined at this stage) can be assessed by comparing the 
separate the — nix mass dependence of m v and m v r obtained in the 
present parametrisation of LaM with the predictions of ChPT. In Fig. 1 
the comparison is done for = and the agreement is fairly good up to 
m K w 800MeV. 

The combination M 2 of f\ and fi 2 , can be split up only by making use 
of the expression of the admixed scalars, therefore the use of one charac- 
teristics of the mixed scalar spectra is unavoidable 22 . 

For the separate determination of /Ltg and fi we have fixed first the mass 
of the fo meson. This meson mass is identified with the heavier eigenvalue 
calculated in the mixing scalar subspace. Since there is the possibility of 
mixing between the normal and the superbound multiplcts we have chosen 
as a second possibility also fixing the trace m 2 + m 2 Q — m 2 ^ + m 2 ^ (see 
Table 1) in the squared mass subspace of admixed scalars. In order to 
test the sensitivity of the results to the scalar masses (which are not very 
accurately known) we have performed the calculations for both alternatives 
with several numerical values. 
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Figure 1. The tree— level kaon mass dependence of m v and m i for = 0. The 
labels refer to the results of ChPT and the predictions of linear sigma model (LctM), 
respectively. 

3. Quasiparticle thermodynamics and phase diagram 

The rcnormalised set of the equations of state and of the sclfconsistent pion 
propagator determine the temperature dependence of the vacuum expecta- 
tion values x and y. The scheme of the perturbation theory applied here 
agrees with the Optimal Perturbation Theory of Chiku and Hatsuda 23 , 
which was renormalised using the approach of 24 . 

The tree level mass of it involves now the thermal mass parameter: 

ml = M 2 (T) + 2(2A + f 2 )x 2 + A.hy 2 + 2gy, (11) 

and all other meson masses to be used in the tadpole integrals below agree 
with the formulas appearing in Table 1 with the replacement — ^ — > 
M 2 (T). If all quantum corrections are condensed into M 2 (T), then the 
tree-level masses of other mesons are expressible through the mass of the 
pion. One might expect that the pion has the lowest mass and therefore 
for M 2 (T) > these squared masses are all positive, which is not the case 
when — Hq < is used in the propagators. We define a physical region of 
x and y where all tree-level mass squares are positive, and thus the one- 
loop contribution of the meson fluctuations to EoS is real. This region is 
most severely restricted by the mass of <x, which strongly decreases near 
the phase transition. We will restrict our attention to the solution of the 
EoS's in the physical region. 
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For the determination of the thermal mass we use the Schwingcr-Dyson 
equation for the inverse pion propagator at zero external momentum. At 
one-loop it receives the contribution H(M(T),p = 0), which is the self- 
energy function of the pion at zero external momentum, plus the coun- 
terterm contribution — ji^ — M (T). We apply the principle of minimal 
sensitivity (PMS) 23 , that is we require that the pion mass be given by its 
tree-level expression: 

n(M(T),p = 0)-/*§-M 2 (T) = 0. (12) 

H(M(T),p) itself is a linear combination of the tadpole and bubble dia- 
grams (the latter not included in the treatment of 12 ), with coefficients 
derived with help of the 4-point and 3-point couplings among mass eigen- 
value fields. 

The self-energy can be represented as a linear combination of tadpole 
integrals, which gives when substituted into Eq. (12): 

a— (7, 7T 

= -M 2 (T)- M 2 + J2 cZJ(m ai (T),T). (13) 

i— 7r, K : rj^rj 

Here c£ are the weights of the renormalised tadpole contributions evaluated 
with different mass eigenstate mesons a, = Cj, 7Tj. The integrals over the 
corresponding propagators are evaluated with effective tree-level masses 
where M 2 (T) replaces — li 2 ,. In this way (13) is actually a gap equation 
which determines the thermal mass parameter, M 2 (T). With help of Eq. 
(11) this equation can be also understood as a gap equation for the pion 
mass (the pion mass is present also in the expressions of I(m ai , T) through 
m ai \): 

a — (T 7 7T 

ml = -fi 2 + 2(2f 1 +f 2 )x 2 +Af 1 y 2 + 2gy + £ J(m ai (T), T) . (14) 

For the determination of the order parameters x, y we solved the two 
renormalised equations of state using the solution of the gap equation in 
the propagator masses: 

-e x - [i\x + 2gxy + 4f lX y 2 + 2(2.h + f 2 )x 3 (15) 
-e y - n 2 y + gx 2 + A.hx 2 y + 4(/i + f 2 )y 3 (16) 
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The quantities t%. and t v a . give the corresponding weights. Ji is the isospin 
multiplicity factor: = 3, Jk — 4, and J v , v i = 1. The coefficients c and 
t were listed in 4 . One finds c„. = J^./x, which ensures that the solution 
for the mass of the pion obeys Goldstone's theorem. 

The solution of Eqs. (14), (15), (16) for given ra„ n%K allows to establish 
the nature of the temperature driven transitions. First order transitions 
are signalled by multivaluedness in the temperature evolution of both the 
non-strange and strange condensates. For large values of the kaon mass, we 
claim that the phase transition is driven by the variation of the non-strange 
condensate, since each of the multiple solutions of the strange condensate 
are very close to each other, and all stay at high values. Subsequent decrease 
of the strange condensate at higher temperature displays only a crossover. 




Figure 2. The temperature dependence in the physical point of: (a) the non-strange 
(x) and strange (y) condensates; (b) the pseudoscalar (8 V ) and scalar mixing angles 
(in the (0-8) basis); (c) the mass of the chiral partners (it, a) and (ao,rj); (d) the mass 
of /o, k, r)' , K mesons. 



In the physical point the system exhibits clear crossover as can be seen 
in Fig. 2. Let us discuss next what happens near the = mx diago- 
nal. In Fig. 3 the gradual deformation of the boundary line is shown when 
four different (m T , m /^-independent conditions (listed on the figure) are 
imposed on the spectra of the admixed scalars. The boundary reaches the 
diagonal in the range m cr it ~ 20 MeV independently of the condition im- 
posed. When moving away from the diagonal for ttlk > the larger is 
the scalar mass scale the farther the boundary goes away from the = 
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Figure 3. Variation of the boundary of first order transitions near and above the m n = 
ma diagonal (note the rather different units on the two axes!) 

axis. Therefore with this method one cannot make any definite statement 
on the location of the tricritical point on the m^-axis defined as the point 
where the transition changes from a discontinuous nature into a continuous 
one. 

In conclusion of this study we find that the assumption of the — tjik 
independence of the mass of the heavier isoscalar-scalar member of the 
nonet, suggested by Gribov's confinement picture leads to a rather unique 
conclusion, that the m n ~ vtlk diagonal crosses the critical curve of the 
chiral phase transitions of the 3- flavor QCD at a rather low value: 15 MeV < 
m C rit < 25 MeV. 
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